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Abstract 

We carried out a study of the properties of the \(f) 4 field solutions. By 
constructing Gaussian wave packets to calculate the S matrix, we show 
that the probability of the vacuum unbroken state transfers to the broken 
state is about 10 -52 . After adding this probability restriction condition 
as modulation factor in the summation of vacuum energy, we thus get a 
result that the vacuum energy density is about 10 -47 GeV 4 , which is exact 
same as the observed dark energy density value, and maybe served as a 
possible explanation of the dark energy. Also our result shows that the 
vacuum energy density is proportional to the square of the universe's age, 
which fits the Dirac large numbers hypothesis. 
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1 Cosmological Constant Problem 

In the current standard model of big bang cosmology, the ACDM model, dark 
energy is a hypothetical form of energy in the space that accelerates the expan- 
sion of the universe. The most recent WMAP [1] [2] observations show that the 
universe is made up of about 74% dark energy, 22% dark matter, and 4% ordi- 
nary matter. The dark energy is measured in the order of about 10 _47 GeV 4 [3J. 
A possible source of the dark energy maybe the cosmological constant. 
The Einstein's gravity field equation is 

1 8ttG 

Rfiv - -^Rg^v + A# MI/ — — j-Tju,, (l) 
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in which A is the cosmological constant, and it can also be written in the form 
of vacuum energy density 



/^vacuumC — ~ • (2) 
07TG 

The general relativity field equation ([1} does not give any physical origin 
of the cosmological constant. We may explain it as the zero-point energy of 
the quantum fields. In quantum field theory, states can be treated as a set of 
harmonic oscillators which poses the zero-point energy 

E = Y,\*»i, (3) 



which is obviously diverge when sum over all states. If we cut off the summation 

G 



at the Planck energy scale Epi an k — J ^fr ~ l0 19 GeV, the zero-point energy 



density will be 

4ir 2 h 3 c 5 

which is more than 10 120 times larger than the measured value of dark energy 
density. To cancel almost, but not exactly, the quantum field theory faces a 
challenge. 

In this article, we try to treat the Higgs field energy as the vacuum energy. 
First we study the static solution of the A0 4 field equation. Then we figure out 
the probability of the vacuum unbroken state transfers to the vacuum broken 
state. Finally with probability modulation factor, we calculate the vacuum 
energy. Our result will be compared to the experimental data. 

2 Solutions Of the A0 4 Field 

The Lagrangian density of the X<p 4 field in the (1+3) dimensions is: 

£ = -i^V-^V-^ 4 , (5) 

in which <\> is the Higgs field, and g„ v = (1,1, 1, 1). The Higgs boson has not been 
found yet in experiments, nor does the Standard Model predict its mass. Varies 
theories and experiments estimate its mass between 115GeV to 180GeV [5]. In 
the Standard Model, with /i 2 < 0, the Higgs boson mass is given by m# = 
\f— 2/i 2 = v2Ai>, where v is the vacuum expectation value of the Higgs field. 



After spontaneous symmetry breaking, v = y — ^- = -J t^jj ~ 246GeV^, fixed 
by the Fermi weak coupling constant Gp . The vacuum unbroken state has 
{4>) = 0, and the vacuum broken state has ((f)) = v. 
From the Lagrangian © , the equations of motion is 

V2( ^^ = ^ + A ^> ( 6 ) 
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which can also be written as 



in which mjj is the Higgs mass, v is the broken vacuum expected value. 

In one dimension, we know equation ^ has static kink/ anti-kink solutions 

cf) = ±v tanh(^-x). (8) 

Because m# ps 115 ~ 180GeV^, slope of (|8j) at a; = is 

dcj)[x) 



d.r 



= 2vm H »Q, (9) 

z=0 



we can replace the 1-dimension kink solution in (jHJ) with a step function in our 
calculation without sacrifice accuracy 

, f v, when x > 0, «2/ \ 2 /, n \ 

^ = { whenx<0, =* A*)=<A (10) 

By study the properties of the Jacobi elliptic functions sn and cn, we can 
get other static solutions of equ. ([6]) in 1-dimension as 

<f>(x) = ±v coth(^-x), (11) 
<j>{x) = ±V2usech(i^|-ar) = iv^u sec(^ra;) = -^j^t, (12) 

4>(x) — ±i\/2 ucsch(z^^-x) = ±V2v csc(^Lx) — r~ 777^- ■ (13) 

v2 v2 sm(-^S-x) 

If 4>(x) — F{Ax) is a solution of the equ. (|6]) in 1-dimension, in which F() 
is a function such as listed above in ©, (ITT1) . (TH21) . (TT3")) and A is a constant, we 
can construct (1 + 3) dimensions solutions of the equ. (JBJ) as 

4>M = F(-^(a ll x li )), (14) 

in which a M are arbitrary constants and automatically sum over \x. For instance, 
it's easy to verify the (1 + 3) dimensions (it,x,y,z) kink-like solution for equ. 
@ as 

<f>{xu) = ±v tanh — - mH = (ia t + a x x + a 2 y + a 3 z) . (15) 

V 2Vflo + a l + fl 2 + fl l / 
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3 Unbroken And Broken States Of the Vacuum 



The Higgs mechanism requires the vacuum unbroken state has {4>) = 0, and the 
vacuum broken state has (4>) — v. Solutions (fT2|) and (|T3|) have domain walls, 
and can only take values above +y/2v or lower than — y/2v, which does not fit 
the (<p) = requirement of the unbroken vacuum state. So we will not use them 
in this calculation. 

To the solutions (JSJ) and (fTTj) , because of their anti-symmetry shape respect 
to = 0, both solutions have the property of (</>) = 0. And because of their 
asymptotic behavior at infinity, if x^ goes like — > oo — > 0, then 4>{x^) goes like 
— > v — » oo. Combination of the solutions ([8]). (fTTj) and transform method (|14| 
can let the 4> field get values from — oo to oo, which satisfied our requirement. 
So we will use static 3-dimensions solution 

<j){x, y, z) = ±v tanh — =JL==(aiX + a 2 y + a 3 z) (16) 
\2y/a{ + aj + a 2 J 

as the initial static state wave function for the unbroken vacuum state. The 
relation <fr 2 ss v 2 is still hold even in the 3-dimension case. That the positive 
or negative sign of the tf> field does not matter, because we will only need </> in 
the following integration. This approximation greatly reduces the calculation 
complexity. 

To distinguish from the unbroken state we will use tp for the spontaneous 
symmetry broken vacuum state. We can construct a static local Gaussian packet 
to represent it, which is centered at oo with a narrow width. Or practically to 
say, it is centered at point b = (b x , b y , b z ), in which (b x , b y , b z ) are the coordinates 
of point b. 

-r2((x-b x ) 2 + ( H -b„) 2 + ( 2 -6 2 ) 2 ) 

ip = ve 5 '- e <P*x+P v y+P*z) ( 17 ) 

in which T 2 is the packet width in momentum space; p x ,p y ,p z are the momen- 
tum of the field; v reflects the role that equ. ((HJ and (fTTj) envelope on the 
Gaussian packet. If values of each of (b x ,b y ,b z ) is very large, (p will be very 
close to v. In fact, since the slope in ([9]) is so large, point b even does not need 
to be very far away from zero. 



4 Packet Width 

We now try to estimate the width of the Gaussian packet in equ. (|T7|) . By 
definition, vacuum is the original point for measuring the energy. So for the 
unbroken vacuum state, energy -E U nbroken = 0. From the Lagrangian©, we 
read out the unbroken vacuum state has m 2 — /i 2 , with /i 2 < 0. From the 
energy-momentum relation 

E 2 =p 2 + m 2 , (18) 
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the momentum will be 

2 

2 i '2 2 2 H l c\\ 

Punbrokcn = ^unbroken — m = ~ M = — ^~ • (.■'■") 

For the vacuum state after spontaneous symmetry breaking, because we have 
moved the original point to (cf>) = v, the energy will be changed. If we think 
momentum is conserved, p U nbroken = Pbroken, and now the broken state has mass 
square as m 2 = m 2 H — — 2/i 2 , then 

7~i2 2 2 3?Tij^- , , 

-^broken = Pbroken + m H = — ^ — • (20) 

From equ. (|18|) . we have _EcLE = pdp. If we treat and c?p approximately 
as and T p for the width of energy and momentum respectively, as explained 
in detail in ref. [3], 

T E E = T p p, (21) 

which gives 

r p = B^il Te = V3r B . (22) 

Pbroken 




Figure 1: List of most significant decay channels branch ratio. (Courtesy from 
CERN website 0). 

Figure (JT|) summarizes the branching fractions of the most important decay 
channels of the Higgs field. The Higgs decay width can be calculated in the 
Standard Model, 



T(H^ff) = 
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Figure 2: Total Higgs width. The dots are simulation results from CMS (Physics 
TDR) for the H — > ZZ* — > 4 leptons final state sub-threshold decay. The solid 
line is the Standard Model prediction. (Courtesy from Tally's talk [6]). 
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which is plotted in Fig. ([2]). When mn < 140GeV, the dominant channels is 
the H — > bb, and the total decay width is in the lOMeV range. 

The results in Fig. ([2]) is arguably for its applicability. When rn# is less 
than 2mw or 2m z, we may research the 'sub-threshold decay', such as H — > 
ZZ*, where one Z boson is on-shell and the second Z boson's mass is off-shell. 
Simulation from CMS (Physics TDR) for the H -> ZZ* ->4 leptons final state 
result is also shown in Fig. © [B] [7]. According to equ. H -> WW 

channel's width should be a little greater than H — > ZZ . But experimentally, 
it's very hard to measure the energy and momentum of the neutrinos in the 
H — > WW — >-2 charged leptons + 2 neutrinos case; compares to the easy be 
detected one in the H — > ZZ — s-4 charged leptons channel. 

For the Higgs mass around ISOGeV, ZZ* channel sub-threshold decay has 
r ~ 2GeV. If we think the WW* channel is about same, then the total width for 



m H = ISOGeV isT E 
for m H = UOGeV. 



AGeV. That means, from ([22]), T p = V3T E » 6.8GeV 



5 Transition Probability Of the Vacuum States 

The lowest order S matrix involves the transition amplitude of the </> 4 interac- 
tion, which made of two incoming vacuum unbroken states 4>n4>i2 in equ. (|16[) . 
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and two out-going broken states tp$\ipfi in equ. (fTTj) . 
4! A 



S if - 



2x2 4 
3Xv 2 



7 (0il(j>i2\<PflVf2) 



if flip f 2 d A x 



3^3/2^4 _ £l±g|+gj ,n<™ 

— — — - — e p x phase x displacement x / dt 

2r„ Jo 



- p 



37r 3 / 2 Aw 4 r 



= 10 



2r3 

-26 



e r p x phase x displacement 



x phase, 



in which the phase and displacement terms are 

phase = e i (p= ( - b ^+ b2 ^+Pv( bl y+ b2 y)+P"( bl '+ b2 ^) ; 



(24) 
(25) 



displacement = exp 



6 2 ,) 2 + (h 



b, 7 \ 



(26) 



The approximate sign in equ. (|24[) means we replace the square of tanh function 
with tanh 2 w 1 as in (fT0|) . So only the Gaussian integrations of ip are left. The 
phase term does not contribute to the following probability calculation. So 
we just omit it. The displacement term reflects the effect that the two out- 
going broken states centered at different points 61 = {bi x ,bi y ,bi z ) and 62 = 
$2x1 bi yi 6 2z ). We just let them to be at same points, so the displacement term 
reaches its maximum value, which is 1. Other values we used are, v = 246GeV, 
A = m 2 H /2v 2 w 1/8, tuh = lSOGeV^ momentum square p 2 — m 2 H /2 is in equ. 
(fl"9|) , T p = \/~3Te — 6.8GeV. The time integral gives out a term, T, which is the 
age of the universe, equals to about 1.3 x 10 10 years = 6.2 x 10 41 GeT/ _1 . Then 
we get the probability that the vacuum spontaneously transfer from unbroken 
states to the broken states, 



probability = =10 



-52 



(27) 



6 Summation Over Possible States With Prob- 
ability 

Vacuum does not like matter, in which matter has so many possible states 
that can oscillate in many frequencies. The vacuum energy can not contain 
contributions from all arbitrary wavelengths, except those physically existed 
states that are permitted by the Hamiltonian, e.g. To calculate the summation 
of states, we made the following assumptions: 

1) The only contribution to the vacuum energy comes from the spontaneous 
symmetry breaking, the Higgs mechanism; which means we do not care any 
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other possible sources. Vacuum has only two states: unbroken state and broken 
state, which will be the only two permitted states that appeared in the sum- 
mation. Because the vacuum unbroken state has the energy E — 0, so the only 
state left in the summation is the vacuum broken state. 

2) The summation needs some modulation factor. For instance, to avoid 
diverge in the thermal radiation expression, in the derivation of the Planck's 
formula of blackbody radiation, we not only add up the energy of photons, but 
also times the Bose-Einstein distribution factor l/(ew — 1). In this paper, we 
chose the probability of the transition from unbroken state to the broken state 
as the suppression factor. 

Finally we sum the energy in phase space and get the vacuum energy density, 

P = 4 ^^3 ~ x Probability = Kr 47 GeF 4 , (28) 

in which the momentum and energy are in equ. (fT9|) and (|20l) . sw (1 to 2) x 
10 2 GeV; dp m T p ; probability term is in equ. (|27p . It is amazing to see that 
our result in (|28|) is exact same as the value that observed in the experiments. 



7 Discussion 



By adding certain conditions in the calculation of the vacuum energy, we get a 
same result to the observed value and maybe served as a possible explanation of 
the dark energy. Higgs has not been found yet in any experiment. The numerical 
estimation of the transition probability is sensitive to the input mass and width 
value. So we need further experimental data to consolidate out calculation. 

When we integrate the S matrix as transition amplitude, we use static wave 
functions. It maybe asked why we have E, Te, p and T p , but only e p x 
and e lpx appeared in the calculation expression, while e~ TBt and e~ tEt did not 
show up? A possible explanation is that E is only occurred after symmetry 
breaking. Its value affects the decay after symmetry breaking, but not related 
to the unbroken state. 

The Dirac large numbers hypothesis [8] [9] [10], made by Paul Dirac in 1937, 
proposed relations of the ratios between some fundamental physical constants. 
The hypothesis has two important consequences: 

1) Matter should be continuously created with time. The mass of the uni- 
verse is proportional to the square of the universe's age: M oc T 2 . According 
to Dirac, the continuous creation of mass maybe either of two models. Additive 
creation, assumes that matter is created uniformly throughout space. Or multi- 
plicative creation, created of matter where matter already exists and proportion 
to the amount already existing. 

2) The gravitational constant, G, is inversely proportional to the age of the 
universe: G oc 1/T. 

In our analysis, the transition amplitude (|24p proportions to T, means that 
the vacuum energy density in equ. ()28[) has the character of p oc T 2 , in which T 
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is the age of the universe. Since the dark energy is about 74% of the universities 
mass, we can roughly to say that equ. (l28l) fits the consequence of the Dirac large 
numbers hypothesis, in case we did not count the expansion of the universe. If 
our calculation is correct, this relation will have a serious effect on the evolution 
of the universe. 

8 Remarks 

The 1-dimension kink/anti-kink solution (|5J| is printed in many books. But we 
had not seen anywhere lists solutions like pT|) . (TT2"j) . (fTB"]) and (TT4")) . For such 
a non- linear equation, we do not know any systematic method to solve it. So 
we just use the dummy way by trying after trying. We note here for further 
suggestion. 
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